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ABSTRACT 

Background: Obsessive-compulsive disorder (OCD) is characterized by maladaptive repetitive behaviors that 
persist despite feedback. Using multimodal neuroimaging, we tested the hypothesis that this behavioral rigidity 
reflects impaired use of behavioral outcomes (here, errors) to adaptively adjust responses. We measured both 
neural responses to errors and adjustments in the subsequent trial to determine whether abnormalities correlate 
with symptom severity. Since error processing depends on communication between the anterior and the poste- 
rior cingulate cortex, we also examined the integrity of the cingulum bundle with diffusion tensor imaging. 
Methods: Participants performed the same antisaccade task during functional MR1 and electroencephalography 
sessions. We measured error-related activation of the anterior cingulate cortex (ACC) and the error-related neg- 
ativity (ERN). We also examined post-error adjustments, indexed by changes in activation of the default network 
in trials surrounding errors. 

Results: OCD patients showed intact error-related ACC activation and ERN, but abnormal adjustments in the post- 
vs. pre-error trial. Relative to controls, who responded to errors by deactivating the default network, OCD pa- 
tients showed increased default network activation including in the rostral ACC (rACC). Greater rACC activation 
in the post-error trial correlated with more severe compulsions. Patients also showed increased fractional anisot- 
ropy (FA) in the white matter underlying rACC. 

Conclusions: Impaired use of behavioral outcomes to adaptively adjust neural responses may contribute to symp- 
toms in OCD. The rACC locus of abnormal adjustment and relations with symptoms suggests difficulty suppress- 
ing emotional responses to aversive, unexpected events (e.g., errors). Increased structural connectivity of this 
paralimbic default network region may contribute to this impairment. 

© 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Obsessive-compulsive disorder (OCD) is characterized by stereo- 
typed and repetitive behaviors that persist despite feedback. These be- 
haviors serve to reduce distress, but are often excessive and not 
realistically connected to the feared outcome they are intended to 
prevent. For example, an individual with OCD may feel compelled to re- 
peatedly check that the door is locked despite having just locked the 
door, 'knowing' that the door is locked, and wanting to stop. In this sce- 
nario, the distressing feeling that the door is not locked persists despite 
clear evidence to the contrary and compels repetitive checking. Such re- 
petitive, maladaptive behaviors may reflect impaired use of feedback 
about outcomes (e.g., a locked door) to adjust emotional (e.g., distress) 
and behavioral (e.g., checking) responses. Here, we used multimodal 
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neuroimaging to test the hypothesis that impaired use of outcomes - 
in this study, errors - to adaptively adjust future responses characterizes 
OCD and contributes to symptoms. Since error processing involves 
(i) recognizing that an error has occurred and (ii) adjusting future re- 
sponses and since deficits in either of these abilities could contribute 
to rigid, repetitive behavior, we examined the neural and behavioral 
markers of each. 

We used an antisaccade paradigm to study error processing. 
Antisaccades require inhibition of the prepotent response of looking to- 
wards a suddenly appearing stimulus and the substitution of a gaze in 
the opposite direction. Antisaccade errors (i.e., looking towards the 
stimulus) reliably elicit neural and behavioral error markers (Agam 
et al., 2011; Belopolsky and Kramer, 2006; Endrass et al., 2007; Klein 
et al., 2007; Nieuwenhuis et al., 2001 ; Polli et al., 2005). We first inves- 
tigated whether individuals with OCD show intact error detection as 
indexed by error self-correction and two extensively characterized 
and highly reliable neural error markers: the error- related negativity 
(ERN) as measured by electroencephalography (EEG) and functional 
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MRI (fMRI) activation of the anterior cingulate cortex (ACC). These error 
markers have been theorized to index error-based reinforcement learn- 
ing (Holroyd and Coles, 2002), mismatch or conflict monitoring (Carter 
and van Veen, 2007; Falkenstein et al., 2000; Yeung et al, 2004), in- 
creased cognitive control (Kerns et al., 2004) and the emotional re- 
sponse to errors (Proudfit et al., 2013). Prior studies of OCD have 
demonstrated an exaggerated ERN (Carrasco et al., 2013; Endrass 
et al., 2008, 2010; Gehring et al., 2000; Johannes et al., 2001; but see 
Nieuwenhuis et al., 2005; Ruchsow et al., 2005; Xiao et al., 2011), in- 
creased rostral ACC (rACC) and/or dorsal ACC (dACC) activation follow- 
ing errors (Fitzgerald et al., 2005, 2010; Maltby et al., 2005; Stern et al., 
201 1 ; Ursu et al., 2003), and inappropriate error signaling on correct tri- 
als in some (Maltby et al., 2005; Ursu et al., 2003), but not all studies 
(Fitzgerald et al., 2005; Gehring et al., 2000). We expected to replicate 
prior findings that exaggerated error signaling correlates with symptom 
severity in OCD (Fitzgerald et al., 2005; Gehring et al., 2000; Ursu et al., 
2003). This would support the theory that inappropriate and exaggerat- 
ed error signaling leads to a pervasive sense of incompleteness and self- 
doubt and triggers compulsions to repeat behaviors ( Maltby et al., 2005 ; 
Pitman, 1987). 

Second, we investigated whether patients with OCD use errors to ad- 
just future responses. To this end we examined post-error slowing ( PES; 
Rabbitt, 1966), or the slowing of responses in trials that follow errors, 
and its neural correlates. Over trials, responses speed up until an error 
is committed (i.e., pre-error speeding; Gehring and Fencsik, 2001; 
Ridderinkhof et al., 2003) and following errors, responses slow down, 
presumably to reduce the probability of another error. (PES does not al- 
ways lead to better performance, however, alternative accounts of its 
function exist (for review see, Danielmeier and Ullsperger, 2011) in- 
cluding that it reflects an orienting response to infrequent events 
(Notebaert et al., 2009).) These adjustments of reaction time (RT) are 
paralleled by changes in activity of the brain's default network (Agam 
et al., 2013). The default network, which is thought to mediate self- 
referential and affective processing, is typically deactivated during ef- 
fortful cognitive performance, presumably reflecting reduced focus on 
the internal milieu (Buckner et al., 2008). During trials immediately pre- 
ceding errors (Agam et al., 2013; Eichele et al., 2008; Li et al., 2007) and 
during error trials (Agam et al., 201 3 ; Polli et al., 2005), there is a relative 
failure of this task-induced deactivation, suggesting that interference 
from internally directed thought culminates in an error. In trials that fol- 
low errors, task-induced deactivation is re-instated (Eichele et al., 
2008), suggesting a shift in focus from the internal milieu back to the 
task-at-hand (Agam et al., 2013). Such dynamic modulations of atten- 
tion and performance in response to outcomes are fundamental to 
adaptive, flexible behavior and we hypothesized that individuals with 
OCD would show abnormal post-error adjustments. We compared RT 
in correct trials immediately preceding an error and correct trials imme- 
diately following an error (i.e., PES), and the corresponding pattern of 
default network activation, as we have done in our prior work (Agam 
et al., 2011, 2013). Use of only pre-error and post-error correct trials 
maximizes sensitivity to changes in performance and activation based 
on error history because it minimizes contamination by global fluctua- 
tions in performance over time that affect comparisons that include all 
correct trials that follow a correct response (Dutilh et al., 2012). We ex- 
pected that abnormal adjustments in the trial after an error in OCD 
would correlate with symptom severity. 

We also examined error positivity or Pe (van Veen and Carter, 2002), 
an event-related potential occurring approximately 300-500 ms fol- 
lowing an error (for review see, Overbeek et al., 2005). The Pe is less 
consistently observed than the ERN and is thought to index error aware- 
ness (Endrass et al., 2007; Nieuwenhuis et al., 2001 ; Wessel et al., 201 1 ) 
and/or the subjective appraisal of errors (van Veen and Carter, 2002) 
and has been associated with PES (Nieuwenhuis et al., 2001 ). 

Finally, we examined the microstructural integrity of the cingulum 
bundle using diffusion tensor imaging (DTI) measures of fractional 
anisotropy (FA). The cingulum bundle contains direct white matter 



connections between the ACC and the posterior cingulate cortex 
(PCC) (Schmahmann et al., 2007), and recent evidence suggests that 
these regions work together to mediate error processing (Agam et al., 
201 1 ). The rACC and PCC are also key anatomical components of the de- 
fault network, whose activation corresponds with PES (Agam et al., 
2013) and may mediate post-error adjustments. We investigated 
whether abnormal cingulate cortex function in OCD is accompanied 
by abnormal structural connectivity. 

In summary, we expected OCD to be characterized by exaggerated 
and inappropriate error signaling and aberrant neural and behavioral 
adjustments following errors. Further, we expected that these abnor- 
malities would predict symptom severity and have an anatomic corre- 
late in the microstructural integrity of the cingulum bundle. 



2. Methods 

2.1. Participants 

Sixty patients from the Obsessive-Compulsive and Related Disor- 
ders Program at Massachusetts General Hospital were enrolled. 
Twenty-seven met the following inclusion criteria and were referred 
for scanning: OCD based on a Structured Clinical Interview for DSM-IV 
(First et al., 1997); Yale-Brown Obsessive Compulsive Scale (Y-BOCS; 
Goodman et al., 1989a,b) total score > 16; no co-morbid Axis I disorder 
with the exception of anxiety disorders and depression; and unmedicat- 
ed or on stable medications for at least 8 weeks. Twenty-one of the 27 
OCD patients referred for scanning completed the study. These partici- 
pants had Y-BOCS scores of 11 ±3 (mean ± SD) obsessions, 12 ± 3 
compulsions, and 23 ± 5 total. Table SI provides co-morbidity and 
medications. 

Twenty healthy control participants, screened to exclude a personal 
history of mental illness (SCID — Non-patient Edition; First et al., 2002) 
and a family history of anxiety disorder, were recruited from the com- 
munity by poster and website advertisements. 

All participants were screened to exclude substance abuse or depen- 
dence within the preceding six months and any independent condition 
that might affect brain function. OCD patients were characterized with 
the Beck Depression Inventory-II (BDI-II, Beck, 1996) and Beck Anxiety 
Index (BAI; Beck etal., 1988; Table SI). The final groups of 21 OCD pa- 
tients and 20 controls did not differ significantly in age, sex, years of ed- 
ucation, handedness based on the modified Edinburgh Handedness 
Inventory (Oldfield, 1971; White and Ashton, 1976), or estimated ver- 
bal IQ based on the Wide Range Achievement Test-Ill Reading portion 
(Wilkinson, 1993) (Table 1 ). The study was approved by the Partners 
Human Research Committee and all participants gave written informed 
consent. 

All participants were included in DTI analyses. To obtain reliable es- 
timates of differences in activation, only participants with a minimum of 
ten usable error trials (Olvet and Hajcak, 2009; Pontifex et al., 2010) 



Table 1 

Means, standard deviations, and group comparisons of demographic data. The Phi value is 
the result of a Fisher's exact test. The z value is the result of a nonparametric Mann- 
Whitney U comparison. 



Subject characteristics Healthy controls OCD t p 

(n = 20) (n = 21) 



Age 


33 ± 11 


33 ± 11 


0.07 


.94 


Sex 


11M/9F 


8M/13F 


<)> = .17 


.35 


Years of education 


17 ± 2 


16 ± 1 


1.33 


.19 


Laterality score (handedness) 3 


65 ± 55 


74 ± 34 


0.62 


.54 


Estimated verbal 10_ b 


113 ± 6 


110 ± 11 


1.10 


.28 



* Laterality scores of — 1 00 and + 1 00 denote exclusive use of left or right hand, respec- 
tively. 

b Two control participants and one OCD patient were not administered the WRAT-III 
because they were non-native English speakers. 
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were included in comparisons of error with correct trials (fMRI: n = 19 
OCD, n = 15 controls; EEC: n = 19 OCD, n = 16 controls), and only par- 
ticipants with a minimum of ten usable pairs of correct pre- and post- 
error trials were included in fMRI analyses of post-error adjustments 
(n = 17 OCD, n = 14 controls). 

2.2. Antisaccade paradigm 

The antisaccade paradigm is described in Fig. 1. The task was pro- 
grammed in MATLAB Psychtoolbox (Mathworks, Natick, MA) and 
consisted of three types of antisaccade trials: hard (40%), easy (50%), 
and fake-hard (10%). Hard trials introduced a distraction during the 
gap, a 3 dB luminance increase of the peripheral squares that mark the 
location of stimulus appearance. Fake-hard trials started with a cue indi- 
cating a hard trial, but were otherwise identical to easy trials (i.e., there 
was no luminance change). They were included as a control condition to 
allow us to examine the effects of a hard vs. easy cue on fMRI activation 
unconfounded by the change in luminance that characterizes hard tri- 
als. As this was not the goal of the present study, error and correct trials 
were combined across trial types for analysis. 

Antisaccade trials were balanced for right and left stimuli. Randomly 
interleaved with the saccadic trials were fixation epochs lasting 2, 4, or 
6 s, which provided a baseline and introduced "temporal jitter" to opti- 
mize the analysis of rapid presentation event-related fMRI data 
(Buckner et al., 1998; Burock and Dale, 2000; Miezin et al., 2000). The 
schedule of events was determined using a technique that optimizes 
the statistical efficiency of event-related designs (Dale, 1999). Each 
task run lasted 5 min 16s and generated an average of 64 antisaccade 
trials and 20 fixation epochs. 

Prior to scanning, participants practiced in a mock MRI scanner, 
were encouraged to respond as quickly and accurately as possible, and 
were told that in addition to the base rate of pay, they would receive 
5c for each correct response. Participants performed six task runs during 
fMRI and eight runs during EEC. The ISCAN fMRI Remote Eye Tracking 
Laboratory (ISCAN, Burlington, MA) recorded eye position during fMRI 
using a 1 20 Hz video camera. During EEC, eye movements were moni- 
tored using two pairs of bipolar EOG electrodes, one vertical (above 




Hard ! ! Hard 



! Fake-hard ! ! ! ! ! 

i i i i i i 

Cue ■ Fixation 1 Gap ■ Stimulus 1 Fixation 1 

0 300 1800 2000 3000 4000 

Time (ms) 

Fig. 1. Antisaccade paradigm. Schematic and timeline of the three conditions: easy, hard, 
and fake-hard. Each trial lasted 4 s and began with an instructional cue (300 ms), either 
a blue or yellow "X" that indicated whether the trial was hard or easy. The mapping of 
cue color to trial type was counterbalanced across participants. The cue was horizontally 
flanked by two white squares of 0.4° width that marked the potential locations of stimulus 
appearance. 10° left and right of center. The squares remained visible for the duration of 
each run. At 300 ms, the instructional cue was replaced by a white fixation ring of 1.3° di- 
ameter at the center of the screen. At 1800 ms, the fixation ring disappeared (200 ms gap). 
At 2000 ms, the fixation ring reappeared at one of the two stimulus locations, right or left 
with equal probability. This was the imperative stimulus to which the participant 
responded by making a saccade in the opposite direction. The ring remained in the periph- 
eral location for 1 000 ms and then returned to the center, where participants were 
instructed to return their gaze for 1 000 ms before the start of the next trial. Fixation epochs 
were simply a continuation of this fixation display. Hard trials were distinguished by a 3 dB 
increase in luminance of the peripheral squares starting during the gap. Except for the 
hard cue, fake-hard trials were identical to easy trials. 



and below the left eye) and one horizontal. Horizontal EOG activity re- 
corded during a brief calibration allowed an estimate of gaze position 
for scoring antisaccades (Endrass et al., 2005, 2007). 

Eye movement data were scored in MATLAB (Mathworks, Natick, 
MA) using a partially automated program. Saccades were identified as 
horizontal eye movements with velocities exceeding 47°/s. The onset 
of a saccade was defined as the point at which the velocity of the eye 
first exceeded 31°/s. Trials with initial saccades in the direction of the 
stimulus were scored as errors. Reaction time (RT) was defined as the 
onset time of the initial saccade relative to the appearance of the stimu- 
lus. Post-error slowing (PES) was defined as the RT difference between 
correct trials immediately following an error trial and correct trials im- 
mediately preceding an error trial. Self-correction latency was defined 
as the latency between the onset of the initial erroneous saccade and 
the onset of the corrective saccade. Error rates were logit-transformed 
before analysis to normalize their distribution. Group differences in 
error rates, rate of error self-correction, RT and PES were assessed 
using between-group t-tests. PES in each group was assessed using a 
one-sample t-test. 

2.3. MRI acquisition 

Images were acquired with a 3 T Siemens Trio whole body high- 
speed imaging device equipped for echo planar imaging (EPI) (Siemens 
Medical Systems, Erlangen, Germany) and a 12-channel head coil. From 
each participant we acquired a structural scan, six functional scans and a 
DTI scan. 

A high-resolution structural scan was acquired in the sagittal plane 
using a 3D RF-spoiled magnetization prepared rapid gradient echo 
(MP-RAGE) sequence (TR/TE/flip = 2530 ms/3.39 ms/7°, FOV = 
256 mm, 176 1.33 x 1 mm in-plane sagittal slices, 1.33 mm thickness). 
Functional images were collected using a gradient echo T2* weighted 
sequence (TR/TE/flip = 2000 ms/30 ms/90°, 32 contiguous horizontal 
slices parallel to the inter-commissural plane, voxel size: 3.1 x 3.1 x 
3.7 mm, interleaved). The functional sequences included prospective 
acquisition correction (PACE) for head motion (Thesen et al., 2000). 
Single-shot echo planar imaging DTI was acquired using a twice 
refocused spin echo sequence with the following parameters: TR/TE = 
7980/84 ms; b = 700 s/mm 2 ; 10 T2 images acquired with b = 0; 60 dif- 
fusion directions; 128 x 128 matrix; 2x2 mm in-plane resolution; 64 
axial oblique (AC-PC) slices; 2 mm (0 mm gap) slice thickness; scan 
duration 9 min, 44 s. 

2.4. pIRI data analysis 

2.4 A. Surface-based projection 

fMRI analyses were conducted on each participant's inflated cortical 
surfaces reconstructed from the MP-RAGE scan using FreeSurfer (http:// 
surfer.nmr.mgh.harvard.edu) segmentation, surface reconstruction, 
and inflation algorithms (Dale et al., 1 999; Fischl et al., 1 999). Functional 
and structural scans were spatially normalized to a template brain 
consisting of the averaged cortical surface of an independent sample 
of 40 adults (Buckner Laboratory, Washington University, St. Louis, 
MO) using FreeSurfer's surface-based spherical coordinate system that 
employs a non-rigid alignment algorithm that explicitly aligns cortical 
folding patterns and is relatively robust to inter-individual differences 
in the gyral and sulcal anatomy of the cingulate cortex. Cortical activa- 
tion was localized using automated surface-based parcellation software 
(Fischl et al., 2004). 

To facilitate comparison with other studies, approximate Talairach 
coordinates were derived by mapping surface-based coordinates back 
to the original structural volume for each of the individuals whose 
brains were used to create the template brain, registering the volumes 
to the Montreal Neurological Institute (MNI305) atlas (Collins et al., 
1994) and averaging the corresponding MNI305 coordinates. These 
coordinates were transformed to standard Talairach space using an 
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algorithm developed by Matthew Brett (http://imaging.mrc-cbu.cam. 
ac.uk/imaging/MniTalairach). 

2.4.2. Functional data analysis 

In addition to on-line motion correction (PACE), functional scans 
were retrospectively corrected for motion using AFNI (Cox and 
Jesmanowicz, 1999), intensity normalized, and smoothed using a 3D 
8 mm FWHM Gaussian kernel. Functional images were aligned to the 
MP-RAGE scan for each participant. Residual head motion for each par- 
ticipant was calculated as the total displacement in millimeters in the x, 
y, and z directions and the total angular rotation in three directions 
(pitch, roll and yaw) from the initial head position for each run as deter- 
mined by the AFNI motion correction algorithm. Displacement and rota- 
tion values were averaged for the six runs of the task and compared 
between groups. 

Finite impulse response (FIR) estimates (Burock and Dale, 2000; 
Miezin et al., 2000) of the event-related hemodynamic responses 
were calculated for each trial type (i.e., pre-error and post-error, error 
and correct) for each participant. This involved using a linear model to 
provide unbiased estimates of the average signal intensity at each 
time point for each trial type without making a priori assumptions 
about the shape of the hemodynamic response. Estimates were comput- 
ed at 12 time points with an interval of 2 s (corresponding to the TR) 
ranging from 4 s prior to the start of a trial to 1 8 s after the start. Tempo- 
ral correlations in the noise were accounted for by prewhitening using a 
global estimate of the residual error autocorrelation function truncated 
at 30 s (Burock and Dale, 2000). 

To test the hypothesis of exaggerated error signaling in OCD, we 
compared activation for correct and error trials at 6 s: the time point 
of maximum error-related activity in a prior antisaccade study (Polli 
et al., 2005). To test the hypothesis of impaired post-error adjustment 
in OCD, we compared activation during the preparatory period of 
post- and pre-error trials at 4 s: the time of peak preparatory ocular 
motor activity in prior studies (Manoach et al., 2007; Polli et al., 2005). 
Group comparisons used t-tests at each vertex. To examine correlations 
between activation and symptom severity in OCD, we regressed activa- 
tion on Y-BOCS scores. 

To correct for multiple comparisons, 10,000 Monte Carlo simulations 
of synthesized white Gaussian noise were run using the smoothing, re- 
sampling, and averaging parameters of the surface-based functional 
analysis to determine the likelihood that a cluster of a certain size 
would be found by chance for a given threshold (p < .05) on the cortical 
surface and to provide cluster-wise probability (CWP) values. 

2.4.3. Region of interest (ROI) definitions 

For error vs. correct comparisons the rACC and dACC were defined 
using automated surface-based parcellation software (Fischl et al., 
2004) that delineated the cingulate cortex and divided it into rACC, 
dACC, and PCC (Desikan et al., 2006). 

For analysis of post-error adjustments a default network ROI was de- 
fined using seed-based functional connectivity analysis of the functional 
data of all participants. A seed in the PCC was defined, based on prior 
studies (Andrews-Hanna et al., 2007; Van Dijk et al., 2010), as a 4 mm 
radius sphere centered at MNI coordinates x = 0, y = — 53, z = 26. 
Pre-processing of the functional data of all participants involved: 1 ) reg- 
istering the motion-corrected fMRI scans to the Montreal Neurological 
Institute (MNI152) atlas using FSL (FMRIB Software Library, www. 
fmrib.ox.ac.uk/fsl); 2) spatial smoothing using a Gaussian kernel of 
6 mm full-width at half maximum; 3) temporal filtering (0.009 Hz to 
0.08 Hz); and 4) removal of spurious or nonspecific sources of variance 
by regression of the following variables: (a) the six movement parame- 
ters computed by rigid body translation and rotation in preprocessing, 
(b) the mean whole brain signal, (c) the mean signal within the lateral 
ventricles, and (d) the mean signal within a deep white matter ROI. The 
first temporal derivatives of these regressors were included in the linear 
model to account for the time-shifted versions of spurious variance. 



Regression of each of these signals was computed simultaneously and 
the residual time course was retained for the correlation analysis. Func- 
tional connectivity maps were created by computing the Pearson corre- 
lation of the signal averaged across the voxels in the PCC seed region 
with that of every other voxel in the brain. Correlation maps for each in- 
dividual were converted to a map of z-scores using Fisher's z transforms. 

Determination of functional connectivity in the averaged group data 
was based on t-tests of the z-scores at each voxel. The default network 
was defined as vertices with positive connectivity that exceeded a 
Bonferroni corrected (for all gray matter voxels based on segmentation) 
threshold of p < 3.3 * 10 -7 , which set the overall probability to .05 
(Table S2). The functional connectivity analysis was conducted in the 
volume and the resulting statistical map was projected onto the cortical 
surfaces (Fig. SI). 

2.5. EEG acquisition and analysis 

EEC was acquired using a 70-channel electrode cap. Electrode im- 
pedances were brought below 20 kfl at the start of each recording ses- 
sion. All signals were identically filtered to 0.1-200 Hz bandpass and 
digitized at 600 Hz. After excluding noisy EEC channels by visual inspec- 
tion of the raw data, EEG data were re-referenced to the grand average. 
Each participant's continuous EEG data were low-pass filtered at 40 Hz. 
Trials with eye blinks were defined by a difference between the maxi- 
mum and minimum voltage of 150 uV or greater at the vertical EOG 
channel and excluded from analysis. 

EEG data, time-locked to saccadic onset, were baseline-corrected by 
subtracting the mean signal during the 100 ms preceding the saccade 
from the 500 ms that followed the saccade. Data for each of two trial 
types (all correct antisaccades and all errors) were averaged for each 
participant. Instead of using data from either FCz of Cz alone, the ERN 
was derived from the average signal across the following 10/20 loca- 
tions: FC1, FCz, FC2, CI, Cz, C2, CP1, CPz, and CP2 for each participant 
to avoid participant exclusion for bad channels. On average, participants 
had 7/9 usable electrodes. Unusable data was seen in five control and 
nine OCD participants at Cz and in five control and two OCD participants 
at FCz. The peak ERN for the combined group of control and OCD partic- 
ipants was identified within the 200 ms following saccadic initiation as 
the point of maximal difference for the error vs. correct waveforms 
(which occurred at 140 ms). The peak ERN for each participant was 
identified as the point of maximal difference within 50 ms on either 
side of the group peak. The Pe was derived using the same set of elec- 
trodes as the ERN (Overbeek et al., 2005) and because a clear peak 
was absent in many participants (particularly those with OCD), Pe was 
defined as the mean difference in the error vs. correct waveforms be- 
tween 300 and 500 ms following the initial saccade. 

2.6. DTI analysis 

Raw diffusion data were corrected for head motion and residual 
eddy current distortion by registering all images to the first acquired 
T2 (b = 0) image, using the FLIRT tool (Jenkinson and Smith, 2001) 
with a 12-df global affine transformation, available through the FSL. 
The diffusion tensor and FA volumes were reconstructed using the stan- 
dard least-squares fit to the log diffusion signal (Basser et al., 1994). FA 
volumes were registered to the high-resolution structural (Tl ) volumes 
for each participant using the T2 (b = 0) volume as an intermediary. 
Inter-participant registration of individual FA maps to the Montreal 
Neurological Institute (MNI305) atlas was performed using each 
participant's Tl structural image, and the resulting transformation was 
applied to individual FA volumes. The MNI-normalized FA volumes 
were smoothed with a 3D Gaussian kernel with 6-mm full-width at 
half-maximum (FWHM). 

We compared FA in OCD and control groups with t-tests at each voxel 
in the cingulum bundle identified using the Jiilich Histological Atlas 
(Eickhoff et al., 2005). Multiple comparisons correction within the 
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cingulum bundle was based on Monte Carlo simulations of synthesized 
white Gaussian noise using the smoothing, resampling, and averaging 
parameters of the DTI analysis to set the overall threshold to p 
< .05. 

3. Results 

3.1. Antisaccade performance 

Table 2 summarizes behavioral results. OCD patients made more er- 
rors than controls (trend in fMRl, significant in EEG). The groups did not 
differ significantly in the proportion of errors committed on hard trials, 
the proportion of self-corrected errors or the latency of correct re- 
sponses. OCD patients showed a trend to have faster corrective saccades 
than controls during EEG. In the EEG data, controls showed significant 
PES, OCD patients showed a trend, and the groups did not differ signif- 
icantly. Neither group showed significant PES in the fMRI data. Given 
the size of PES (13 ms on average for controls measured by EOG), the 
faster sampling rate of EOG (600 Hz) compared to the fMRI video cam- 
era (120 Hz) may have been necessary to detect PES. In OCD, neither 
error rate nor PES correlated with Y-BOCS scores. 

3.2. Error markers 

OCD and control groups did not differ in residual head motion during 
fMRI (displacement: control 1.9 ± 0.8 mm, OCD: 2.0 ± 0.5 mm, t(37) = 
0.28, p = .78; rotation: control 0.8 ± 0.5°, OCD 0.7 ± 0.35°, t(37) =1.02, 
p = .31 ). Both groups showed significant error-related dACC activation 
and did not differ significantly (Fig. 2A, Table 3A). Similarly, both groups 
showed a significant ERN (control: t(15) = 7.73, p < 10~ 5 ; OCD: 
t(18) = 6.10, p < 10~ 5 ) that did not significantly differ in peak amplitude 
(control: 3.59 ± 1.85 uV; OCD: 3.89 ± 2.77 uV; t(33) = 0.43, p = .67; 
Fig. 2B) or latency (control: 137 ± 29 ms; OCD: 126 ± 28 ms; t(33) = 
1.12, p = .27). The ERN was derived using the average signal across 
nine electrodes. Considering the data from either Cz or FCz alone or 
their average, the ERN did not differ significantly between groups and 
was qualitatively similar to the results averaged over all nine electrodes. 
In the averaged group data, the ERN peaked -140 ms, which is within 
the range of other ERN studies that measured the onset of movement 
with EMG or saccadic measurements (e.g., Gehring and Fencsik, 2001; 
Gehring et al., 1993; Nieuwenhuis et al., 2001), rather than its outcome 
(e.g., button presses), which results in shorter ERN latency (see Agam 
et al., 201 1 for review; Carter and van Veen, 2007; Gehring et al., 2012). 
In both groups, the ERN was maximal at electrode CPz (Fig. S2). Finally, 
the Pe showed a trend to be reduced in OCD patients relative to control 



participants (F(l,33) = 3.98, p = .06; Fig. 2). There were no significant 
correlations of error-related activation in the ACC, ERN peak amplitude, 
or Pe with Y-BOCS scores in OCD. 

3.3. Post-error adjustments 

Control participants showed the expected pattern of default net- 
work activation: a failure of task-induced deactivation in the pre-error 
trial and a reinstatement of significant task-induced deactivation fol- 
lowing an error (Fig. 3A, Table 3B). OCD participants showed an oppo- 
site pattern: significant task-induced deactivation in the pre-error trial 
and a relative failure of task-induced deactivation after the error. This 
pattern of post-error adjustment (post- vs. pre-error trials) significantly 
differed between groups in the left rACC and PCC nodes of the default 
network. OCD patients also showed significantly greater activation of 
the left PCC than controls in the post-error trial vs. fixation. Greater ac- 
tivation of the right rACC in post-error trials in OCD correlated with 
more severe compulsions and total Y-BOCS scores (Fig. 3B; Table 3B). 

3.4. DTI 

Relative to control participants, OCD patients showed increased FA 
in the cingulum bundle (Fig. 4). The peak of the group difference 
corresponded to the right rACC region in which post-error activation 
correlated with symptoms (MNI coordinates: x = 6, y = 30, z = — 4; 
clusterwise probability (CWP) = 0.001). No other regions showed 
group differences that surpassed correction for multiple comparisons 
in the entire brain. FA in the cingulum bundle did not correlate with 
Y-BOCS scores, error-related dACC activation, ERN, or post-error 
adjustment. 

3.5. Control analyses 

3.5.1. Effects of EEG saccadic artifacts on the ERN 

As in prior antisaccade studies (Endrass et al., 2005, 2007; 
Nieuwenhuis et al., 2001 ), the horizontal EOG for saccades to the left 
and right were almost identical in amplitude but of opposite polarity 
(Fig. S3A). Because there were an equal number of stimuli to the right 
and left and the proportion of errors to the right and left did not differ 
within individuals (control: t(15) = 1.09, p = .29; OCD: t(18) = 0.51, 
p = .62), the effects of the horizontal component of the eye movements 
are effectively eliminated in the averaged correct and error waveforms. 
The remaining vertical component of the saccadic artifact peaks at 
19 ms and has the same polarity for left and right (Fig. 2B). It con- 
sists of a positive deflection that returns to baseline -100 ms before 



Table 2 

Means, standard deviations, and group comparisons of task performance. 







Healthy controls 


OCD 




t 


P 


Error rate 3 


fMRI 


12 ± 11% 


14 ± 


9% 


1.67 


.10 




EEG 


13 ± 12% 


16 ± 


8% 


2.10 


.04 


Percent of errors committed on hard trials 3 


fMRI 


72 ± 16% 


70 ± 


15% 


0.36 


.72 




EEG 


75 ± 17% 


70 ± 


13% 


1.54 


.13 


Rate of error self-correction 


fMRI" 


99 ± 1% 


98 ± 


1% 


0.57 


.57 




EEG C 


99 ± 0% 


98 ± 


1% 


0.90 


.37 


Correct response latency 3 


fMRI 


283 ± 11 ms 


267 ± 


10 ms 


1.05 


.30 




EEG 


258 ± 10 ms 


237 ± 


11 


1.35 


.18 


Self-correction latency 


(MRI b 


179 ± 9 ms 


180 ± 


8 ms 


0.04 


.97 




EEG C 


194 ± 8 ms 


174 ± 


8 ms 


1.79 


.08 


Post-error slowing 


fMRI" 


5 ± 19 ms 


6 ± 


27 ms 


0.12 


.91 


(post-error — pre-error latency) 


EEG b 


13 ± 19* 


6± 


13 ms 


1.42 


.17 



Based on 20 control and 21 OCD participants (i.e., the entire sample prior to exclusion based on the number of usable error trials). 
Based on 15 control and 19 OCD participants. 
Based on 18 control and 19 OCD participants. 
Based on 14 control and 17 OCD participants. 

Within-group paired comparison of post and pre-error latency: t(14) = 2.67, p = .02. 
Within-group paired comparison of post and pre-error latency: t(18) = 1.85, p = .08. 
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Fig. 2. Neural responses to error commission. A. Pseudocolor statistical maps of error-related fMRl activation at 6 s in the error vs. correct contrast are displayed on inflated medial cortical 
surfaces. The rACC and dACC ROIs are outlined in black. Gray masks cover subcortical regions in which activation is displaced in surface-based analyses. In the first two rows, warm colors 
indicate greater activation on error than correct trials. The third row shows the group comparison of error-related activation. The middle columns show hemodynamic response functions 
for each condition (correct: black; error: red) and for the error-correct difference for each group (Control: solid line; OCD: dashed line) for the left and right dACC ROIs. B. Plots in the first 
two rows show grand average EEG waveforms, with standard error lines, for correct (black) and error (red) trials, time-locked to the onset of the saccade (0 s), for control (first row) and 
OCD (second row) participants. The first peak after saccadic onset is eye movement artifact, which subtracts out of the difference waveform for error and correct trials, shown in the third 
row for each group. The arrow denotes the approximate time of peak ERN, which did not differ significantly by group. 



the peak of the ERN and is mostly eliminated when the correct and 
error waveforms are subtracted to derive the ERN. It therefore 
does not contribute to the ERN. Nor does artifact from corrective sac- 
cades on error trials substantially affect the ERN. At the same EEG lo- 
cations used to measure the ERN, we quantified the effect of the 
corrective saccade (Fig. S3B). The slope of the resulting signal, 
time-locked to saccadic onset, is smaller that that of the ERN and is 
similar for OCD and control groups. In addition, corrective saccades 
are jittered in time relative to the error, which reduces their effect. 
Finally, on average, corrective saccades are initiated later than the 
peak ERN (controls: 54 ± 33 ms; OCD: 34 ± 35 ms Fig. S3C). 

3.5.2. Comparison of the correct EEC waveform in OCD vs. controls 

A two-way ANOVA with factors trial type (correct, error) and group 
(control, OCD) on the average potentials between 90 and 190 ms after 
the response (i.e., the time window used for defining individual ERN 
peaks), using the same set of nine electrodes as the ERN analysis, 
showed the expected strong main effect of trial type (F(l,33) = 22.05, 
p< 10 -4 ), but no effect of group (F(l,33) = 0.33, p = .57) and no inter- 
action of group with trial type (F(l,33) = 0.20, p = .65). Neither cor- 
rect (t(33) = 0.17, p = .87) nor error (t(33) = 0.76, p = .45) 
waveforms differed between groups indicating that our non- 
replication of the finding of an exaggerated ERN in OCD (e.g., Gehring 
et al., 2000) was not due to an abnormally high amplitude correct 
waveform. 

3.5.3. Effects of group differences in error rate on the ERN, dACC activation 
and neural indices of post-error adjustment 

Error rates often inversely correlate with ERN and dACC activa- 
tion (Holroyd and Coles, 2002; Polli et al, 2008; Santesso et al., 
2006). Regression of activation (or ERN) on error rate with a group 
by error rate interaction showed that this was true in the present 



study for left and right dACC activation (left: r(32) = -0.39, p = 
.02; right: r(32) = — 0.47; p = .005) and that there were no group 
differences in the slopes of these relations (p's > .75). Error rate was 
not significantly correlated with the ERN (r(33) = —0.10, p = .55) 
and there was no group difference in slopes (p = .71 ). While OCD par- 
ticipants made more errors than controls, this difference was eliminat- 
ed after participants were excluded from the functional analyses based 
on having made an insufficient number of errors (fMRI: control 15 ± 
11%, OCD 14 ± 9%, t(32) = 0.11, p = .91; EEG: control 16 ± 11%, 
OCD 16 ± 8%, t(33) = 0.79, p = .44) and therefore does not account 
for the non-replication of findings of exaggerated error markers in 
OCD. A vertex-wise regression of post-error adjustment (contrast 
between post- vs. pre-error trials) on error rate with a group by error 
rate interaction term did not reveal significant relations within the 
default network ROl or the entire cortical surface and this did not differ 
by group. 

3.5.4. Contributions of medication, depression, and anxiety to ERN and 
dACC activation 

Seven OCD participants took various antidepressant medications, 
and one also took memantine (Table SI). Medicated patients did not 
significantly differ from medication-free patients in error-related dACC 
activation (left: t(17) = 0.72, p = .48; right: t(17) = 1.03, p = .32). 
Nor did post- vs. pre-error adjustments in default network activation 
differ significantly between medicated and unmedicated patients 
(t(15) = 0.41, p = .69; Fig. S4A). ERN amplitude was higher in 
medicated patients (Fig. S4B), but this did not reach significance 
(t(17) = 1.68, p = .11). 

BDI-II and BAI symptom rating scale scores in OCD participants did 
not significantly correlate with either error-related dACC activation or 
error-related changes in the default network (i.e., there were no sig- 
nificant clusters). Nor did the ERN correlate with these ratings (BDI-II: 
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Table 3 

Maxima and locations of significant clusters of activation. p-Values ( — base — 10 log) are provided for the most significant vertex. A. Error-related flvlRI activation in the ACC ROls. 
Clusterwise probabilities (CWPs) are based on correction within the entire ACC. B. Default network activation in trials surrounding errors and correlations with Y-BOCS scores. CWPs 
are based on correction within the default network ROI. BA = Brodmann Area. 



Cortical region 



Size (mm ) 



Direction of effect 



Talairach coordinates 



BA 



Max. p-value 



CWP 



Error vs. correct 
Control (n = 15) 
Right dACC 
Right dACC 
OCD (n = 19) 
Left dACC 
Right dACC 
Control vs. OCD 
None 



272 
179 



1068 
969 



Error > correct 
Error > correct 



Error > correct 
Error > correct 



13 
5 



17 
15 



23 
19 



31 
22 



22 
28 



32 
24 



32 
32 



2.45 
2.11 



4.84 
4.16 



.0001 
.03 



.0001 
.0001 



Pre-error vs. fixation 

Control (n = 14) 
None 

OCD (n = 17) 
Left rACC 2365 
Left PCC 784 
Left angular gyrus 551 
Right rACC 2661 
Right PCC 1046 
Right middle temporal gyrus 703 
Right angular gyrus 578 

Control vs. OCD 
None 

Post-error vs. fixation 
Control 

Left rACC 1331 
Left PCC 912 
Left angular gyrus 518 
Right rACC 1757 
Right PCC 577 
OCD 

Right middle temporal gyrus 521 
Control vs. OCD 

Left PCC 584 
OCD, relations with Y-BOCS scores 

Total 

Right rACC 475 
Obsessions 

Right rACC 300 
Compulsions 
Right rACC 558 
Post-error vs. pre-error (adjustment) 
Control 

Left rACC 333 
Left rACC 289 

OCD 
None 

Control vs. OCD 
Left rACC 950 
Left PCC 358 



Negative 
Negative 
Negative 
Negative 
Negative 
Negative 
Negative 



Negative 
Negative 
Negative 
Negative 
Negative 

Negative 

OCD > control 

r> 0 
r> 0 
r> 0 



Negative 
Negative 



OCD > control 
OCD > control 



-2 
-7 
-48 
4 
5 
60 
46 



-10 
-15 
-40 
10 
6 

52 

-16 



11 
10 
14 



-11 
-14 



-12 
-12 



31 

-38 
-57 
31 
-31 

-34 
-57 



41 

-43 
-64 
46 

-58 



-44 

36 
35 
39 



37 
43 



39 
-51 



0 

32 
41 
3 
32 
-5 
42 



20 
36 
39 
12 
34 

-25 

33 



-7 
28 



24 
31 
39 
24 
31 
21 
39 



32 
31 
39 
9 
31 

38 

31 



32 
32 
32 



32 
32 



32 
31 



3.37 
3.23 
2.17 
3.54 
3.26 
3.11 
2.48 



3.05 
2.58 
3.27 
3.75 
2.94 

2.50 

2.08 



2.74 
2.58 
2.56 



3.41 
2.15 



2.78 
2.25 



.0001 

.0001 

.003 

.0001 

.0001 

.001 

.005 



.0001 
.0001 
.004 
.0001 
.005 

.01 

.002 



.02 
.16 
.008 



.008 
.02 



.0001 
.04 



Not significant, but shown in Fig. 3. 



r(17) = -0.27, p = .26; BAI: r(17) = -0.02,p = .94). These findings 
suggest that neither medication nor the presence of anxious or depres- 
sive symptoms substantially affected our findings. 

4. Discussion 

OCD participants showed aberrant use of errors to adaptively adjust 
neural responses despite intact neural and behavioral indices of error 
detection (the ERN, error-related dACC activation and the rate of error 
self-correction). Control participants showed the expected pattern of 
adjustments based on prior work: a failure to deactivate the default net- 
work in pre-error and error trials (Agam et al., 2013; Eichele et al., 2008; 
Li et al., 2007; Polli et al., 2005) and a reinstatement of task-induced de- 
activation in post-error trials (Eichele et al., 2008). This pattern suggests 



that the error and pre-error trial are characterized by interference from 
self-referential thought and that the error prompts a shift in focus back 
to the task-at-hand in the subsequent trial. In contrast, OCD patients 
showed the opposite pattern, increased default network activation in 
the post-error compared to the pre-error trial and differed significantly 
from controls in this regard. This pattern suggests that rather than 
responding to errors by turning attention away from the internal milieu 
and back to the task-at-hand (Agam et al., 2013), the error prompted 
greater self-referential thought. This finding is consistent with previous 
reports that OCD participants show increased default network activa- 
tion, specifically in the ventromedial prefrontal cortex, on error trials 
(Stern et al., 201 1 ) and increased functional connectivity of the default 
network with the fronto-parietal attention network (Stern et al., 
2012). These findings support the hypothesis that OCD is characterized 
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Y-BOCS Compulsions 

Fig. 3. Preparatory fMRI activation of the default network and its relation to symptom se- 
verity. Pseudocolor statistical maps of activation at 4 s are displayed on inflated medial 
cortical surfaces. The default network ROI is outlined in black. Gray masks cover subcorti- 
cal regions in which activation is displaced in surface-based analyses. A: Preparatory fMRI 
activation in the group comparison of the post-error vs. pre-error contrast at 4 s. Blue in- 
dicates regions showing greater activation for OCD than controls in the post- vs. pre-error 
trial. The left graph shows the differences in hemodynamic responses for post-error minus 
pre-error trials for each group at the vertex with the maximum group difference in the 
post- vs. pre-error contrast. The right graph shows the means and standard errors for pre- 
paratory activation at this same vertex in the pre-error trial, error trial, and post-error trial. 
B: Relations between activation on post-error trials and Y-BOCS total, obsession, and com- 
pulsion scores. The scatter plots correspond to the maximum vertices in the correlation. In 
accord with recent practice in fMRI research, r-values from the maximum vertices are not 
reported since these are based on the present data and may therefore inflate the true cor- 
relation by adding the effect of random variability (e.g., Kriegeskorte et al., 2009; Vul and 
Pashler, 2012; Vul et al, 2009). The correlation with obsessions did not meet the cluster- 
wise probability threshold for significance. 



by difficulty disengaging from evaluative processes when errors occur, 
which interferes with redirecting attention to the task-at-hand. The 
present finding of increased FA of rACC white matter in OCD, suggestive 
of increased default network structural connectivity, complements the 
findings of increased default network functional connectivity, although 
the target regions of the presumed increased rACC structural connectiv- 
ity are unknown. Finally, the correlation of greater post-error activation 
of the rACC node of the default network with the severity of compulsions 




Fig. 4. DTI results. Pseudocolor statistical map of group differences in fractional anisotropy 
(FA) displayed on the MNI152 template brain. Red indicates higher FA in controls and blue 
indicates higher FA in OCD. Crosshairs denote the voxel of maximal significance. 



in OCD suggests that maladaptive neural responses to behavioral out- 
comes contribute to symptoms. 

The severity of compulsions in OCD correlated with increased rACC 
activation in the post-error trial, but not with the post- vs. pre-error ad- 
justment as predicted. Thus, this result should be considered prelimi- 
nary. This correlation raises the possibility that more severely affected 
patients are simply more internally focused throughout the task and 
not just after an error. This is unlikely since default network activation 
in neither pre-error trials nor averaged across all correct trials correlat- 
ed with symptoms, even without the cluster-wise correction for multi- 
ple comparisons. This relation may reflect that more symptomatic 
individuals dwell more on aversive and unexpected outcomes such as 
errors, or that aberrant neural responses to behavioral outcomes con- 
tribute to compulsive behavior. Alternatively, both aberrant adjust- 
ments to errors and compulsion severity may be influenced by a third 
factor, such as increased structural connectivity of the rACC as indexed 
by increased FA in the underlying white matter, which could either be 
a trait or reflect experience-dependent plasticity (Scholz et al., 2009). 
The cross-sectional correlative data of the present study does not 
allow us to resolve these important issues. 

We also expected aberrant neural adjustments to errors to lead to 
abnormal behavioral adjustments and poorer performance. Neural ad- 
justments did not correlate with PES or error rate. In the EEG latency 
data, controls showed significant PES, OCD patients showed a trend, 
and although the magnitude of PES in controls was more than twice 
that of OCD patients, the group difference did not reach significance. 
PES has been shown to be relatively small or absent on tasks with longer 
stimulus-to-response intervals as is the case in the present study 
(Danielmeier and Ullsperger, 2011; Jentzsch and Dudschig, 2009). 
Prior studies of OCD are inconsistent showing increased (Fitzgerald 
et al., 2005), intact (Endrass et al., 2008, 2010; Grundler et al., 2009; 
Hajcak and Simons, 2002; Santesso et al., 2006), and reduced (Liu 
et al., 2012) PES. With regard to antisaccade performance, considering 
the entire sample (i.e., before exclusion based on error-rate) OCD pa- 
tients had a higher error rate than controls, consistent with some 
(Lennertz et al., 2012; Rosenberg et al., 1997; Tien et al., 1992) but not 
all (Jaafari et al., 2011 ; Spengler et al., 2006) studies and with other ev- 
idence of impaired response inhibition in OCD (Chamberlain et al., 
2005; Menzies et al., 2007). Neural responses are likely to be more 
sensitive and specific indices of adjustments to errors than behavior. 

Our findings of intact ERN and error-related dACC activation in OCD 
are discrepant with some prior work showing exaggerated error re- 
sponses. Recent meta-analyses (Mathews et al., 2012; Moser et al., 
2013) suggest that whether an exaggerated ERN is found depends on 
the task used and the type of symptoms present. Many prior studies 
used flanker tasks, and we cannot exclude the possibility that the dis- 
crepancies seen in the present study stem from the different movement 
modality or task requirements of antisaccades or that in both groups 
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virtually all errors were self-corrected. Reviews of the literature (Agam 
et al., 2011; Holroyd and Coles, 2002), however, indicate that the ERN 
and dACC activation are 'generic' indices of error processing that are 
present for errors in a range of contexts including errors committed 
with the feet (Holroyd et al., 1998), eyes, and hands in a variety of 
tasks. A prior study that compared ERNs generated by button presses 
and saccades found them to be similar in morphology, amplitude and 
scalp topography (Van 't Ent and Apkarian, 1999). In addition, in our 
prior antisaccade study, the ERN was compatible with those described 
in the literature with regard to scalp distribution, latency, and morphol- 
ogy (Agam et al., 201 1 ). Thus, it is unlikely that antisaccades produce a 
fundamentally different ERN than other tasks used to study error pro- 
cessing. The present findings are consistent with the meta-analyses in 
suggesting that findings of exaggerated error markers in OCD may be 
task-dependent (Mathews et al., 2012; Moser et al., 2013). A caveat to 
this explanation of the discrepancies in our finding of a normal ERN is 
that in the meta-analysis of Matthews and colleagues, increased ERNs 
in OCD were seen for response conflict tasks, of which the antisaccade 
is an example, and, consistent with recent study (Riesel et al., 2014), 
are also seen across all symptom dimensions, except hoarding, which 
was not present in our sample. 

A related concern is that the discrepancies may arise as a function of 
a deficit in volitional saccade generation in OCD. OCD patients and their 
syndromally-unaffected first-degree relatives have been shown to have 
longer latencies of volitional saccades (Kloft et al., 2013). In the present 
study, the latency to generate correct volitional antisaccades was actual- 
ly non-significantly shorter for OCD than control participants suggesting 
intact generation of volitional saccades, despite the deficit in response 
inhibition (i.e., increased antisaccade error rate in the OCD sample 
prior to exclusion). Given that the latency for correct antisaccades and 
the neural markers of antisaccade errors were intact in OCD, it is unlike- 
ly that a deficit in generating volitional saccades accounts for our find- 
ings of intact ERN and error-related dACC activation in the context of a 
specific impairment in post-error neural adjustments. 

Another potential source of inconsistency with the literature is that 
the ERN is variably defined based solely on error trials or on the differ- 
ence between error and correct trials. If the waveform for correct trials 
is also more negative, as is sometimes seen in OCD (Endrass et al., 
2008, 2010; Hajcak and Simons, 2002), difference waveforms are less 
likely to differentiate between groups. This was not the case in the pres- 
ent study since neither the correct nor the error waveforms differed in 
OCD. In the present study, we derived the ERN using the average signal 
across nine electrodes instead of using data from either FCz or Cz alone 
to avoid excluding participants based on bad channels. While we cannot 
exclude the possibility that missing data from either FCz or Cz in some 
participants affected our results, when the data from either Cz or FCz 
was considered alone or when their average was considered, the ERN 
did not differ significantly between groups and was qualitatively similar 
to the results from the data averaged over nine electrodes. 

Methodological differences may have also contributed to the appar- 
ent discrepancy with prior findings of enhanced error-related ACC acti- 
vation in OCD (Fitzgerald et al., 2005, 2010; Maltby et al., 2005; Stern 
et al., 2011; Ursu et al., 2003). Prior studies used analysis techniques 
that assume a shape to the hemodynamic response. A single assumed 
model is unlikely to be valid across all brain regions and stimulus 
types (Duann et al., 2002) and inaccurate models may misattribute 
activation to adjacent events (Manoach et al., 2003). For this reason, it 
is possible that increased ACC activation on error trials in prior studies 
reflects greater activation of the default network prior to error commis- 
sion (Polli et al., 2005) rather than an exaggerated response to the error. 
In the present study, the relatively long cue-to-target interval of each 
trial and the use of FIR models, which make no a priori assumptions 
about the shape of the hemodynamic response, allow more accurate 
modeling of preparatory vs. error-related activation. 

Other potential sources of discrepant findings include medication 
use in the OCD patients and sample size. In the present study, the subset 



of medicated OCD patients did not differ in dACC activation, but showed 
a non-significantly higher ERN, as was also seen in an earlier study that 
showed a normal ERN in OCD (Nieuwenhuis et al., 2005). Medications 
are unlikely to account for previous findings of exaggerated ERNs in 
OCD, however, since this has been reported in OCD patients regardless 
of medication use (Endrass et al., 2008; Stern et al., 2010) and in 
syndromally-unaffected first-degree relatives of OCD patients (Riesel 
et al., 201 1 ). Nor is the non-replication likely to reflect a lack of power 
due to the relatively small sample size. The between-group effect sizes 
of the present study are so small (0.24 in the left dACC, 0.10 in the 
right dACC, and 0.09 for the ERN) that samples of well over 100 partic- 
ipants per group would be required to detect group differences with 
80% power. 

The relatively small sample size may have limited our power to eval- 
uate the Pe, which is less well-characterized than the ERN and less con- 
sistently observed. The Pe showed a trend to be reduced in the present 
study of OCD, but not in prior work (Endrass et al., 2008, 201 0; Ruchsow 
et al., 2005). Unlike the ERN, which is present regardless of whether er- 
rors are perceived, the Pe is present only for perceived errors and is 
thought to index error awareness (Endrass et al., 2007; Nieuwenhuis 
et al., 2001 ). While the literature lacks a clear consensus about its func- 
tion (Overbeek et al., 2005), the Pe has also been associated with short- 
term performance adjustments such as PES (Nieuwenhuis et al., 2001 ). 
Given that in the present study the Pe reduction was not expected, did 
not reach statistical significance and did not correlate with PES, we are 
reluctant to speculate about how an abnormal Pe in OCD may relate to 
our main findings of aberrant neural adjustments to errors in the default 
network. 

Finally, we expected to replicate correlations of increased error signal- 
ing with symptom severity in OCD (Fitzgerald et al., 2005; Gehring et al., 
2000; Ursu et al., 2003). Neither ERN amplitude nor error-related dACC 
activation correlated with symptoms. These null findings are consistent 
with a prior study that reported no relations of ERN with symptom sever- 
ity in OCD (Riesel et al., 201 1 ). This study also found an increased ERN in 
syndromally-unaffected first-degree relatives leading the authors to sug- 
gest that overactive error monitoring is an endophenotype of OCD that is 
independent of the presence of clinical symptoms. In the present study, 
we instead found that increased rACC activation in trials that followed er- 
rors correlated with the severity of compulsions. As this is the first report 
of abnormal post-error neural adjustments, there is no evidence bearing 
on whether this represents a state or trait marker. 

In summary, the present findings demonstrate aberrant neural ad- 
justments to errors in OCD. This suggests an impairment in adaptively 
and dynamically responding to outcomes. Moreover, greater default 
network activation following an error correlated with the severity of 
compulsions. The rACC locus of the aberrant adjustment and relations 
with symptoms suggests difficulty suppressing emotional evaluative 
responses to aversive and unexpected outcomes such as errors. This 
may preclude more adaptive responses that would redirect attention 
to the task-at-hand and minimize the probability of error recurrence. 
Increased FA in rACC white matter suggests increased structural 
connectivity of this paralimbic default network region as an anatomic 
correlate of abnormal adjustment. 

Supplementary data to this article can be found online at http://dx. 
doi.org/1 0.1 01 6/j.nicl.201 4.06.002. 
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